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ABSTRACT 



We have analyzed the properties of the Ha and [Nii]/16583 rest-frame optical emission lines of a sample of 53 intensely star forming 
galaxies at z=1.3 to 2.7 observed with SINFONI on the ESO-VUT. Similar to previous authors, we find large velocity dispersions 
in the lines, cr=few 10-250 km s"'. Our data agree well with simulations where we applied beam-smearing and assumed a scaling 
relation of the form: velocity dispersion is proportional to the square root of the star-formation intensity (star-formation rate per unit 
surface area). We conclude that the dispersions are primarily driven by star formation. To explain the high surface brightness and 
optical line ratios, high thermal pressures in the warm ionized medium, WIM, are required (P/k> 10^-10' K cm"'). Such thermal 
pressures in the WIM are similar to those observed in nearby starburst galaxies, but occur over much larger physical scales. Moreover, 
the relatively low ionization parameters necessary to fit the high surface brightnesses and optical line ratios suggest that the gas is not 
only directly associated with regions of star formation, but is wide spread throughout the general interstellar medium. Thus the optical 
emission line gas is a tracer of the large scale dynamics of the bulk of the ISM. 

We present a simple model for the energy input from young stars in an accreting galaxy, to argue that the intense star-formation is 
supporting high turbulent pressure, which roughly balances the gravitational pressure and thus enables distant gas accreting disks to 
maintain a Toomre disk instability parameter Q~l. For a star formation efficiency of 3%, only 5-15% of the mechanical energy from 
young stars that is deposited in the ISM is needed to support the level of turbulence required for maintaining this balance. Since this 
balance is maintained by energy injected into the ISM by the young stars themselves, this suggests that star formation in high redshift 
galaxies is self -regulating. 
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1. Introduction 



galaxies: formation — galaxies: kinematics and dynamics — galaxies: ISM — galaxies: star 



Galaxies exhibit a wide range of phenomena, some of which 
appear to have been more important in the distant Universe 
than today. The global co-moving star-formati o n rate was 
higher by an order of magnitude (iMadau et al. 19961 ), and 
so were the spec ific star-formation rates (Elbaz et al.l l2007t 
iDaddi et aP l2007h . Morphologies were increasingly irregular, 
an observatio n that cannot be explained with extinction alone 
dConselice et al. 2008) . Processes like starburst-driven outflows 
are o bvious and ubiquitous in samples of hig h redshift galaxies 
(e.g., IShaplev et aLll2003l: ISteidel et alj|2010l) compared to onl y 



a small fraction of galaxies today jLehnert & Hecknianl 1996h. 
The m erger rate is also likely to be much higher (Ide Ravel et al.l 
|2009|) . The gas infall rates of nearby galaxies are apparently 
low, but they were likely much higher in the early Universe 
dCroton et al.ll2006h . Today, the Universe can be considered in 
the "age of secular evolution", but until about 8 Gyrs ago it was 
in the "age of feedback and self-regulation". It is the mix of 
physical processes during their early evolution that shaped the 
characteristics of the ensemble of the galaxies we observe today. 



* Data obtained as part of program IDs: 074.A-9011, 075.A- 
0318, 075.A-0466, 076.A-0464, 076.A-0527, 076.B-0259, 077.B-0079, 
077.B-0511, 078.A-0055, 078.A-0600, 079.A-0341 and 079.B-0430 at 
the ESO-VLT. 



Understanding this complexity and the relative contribution 
of the physical processes that shaped distant galaxies is required 
before we can consider our understanding of galaxies in any 
sense complete. With the higher activity levels of galaxies at 
z~l-5 relative to today it is difficult to gauge the importance 
of each process by mere analogy. This requires careful, direct 
observations of distant galaxies. 

Of significant concern in this pursuit is that obtaining de- 
tailed spatial information from distant galaxies is hampered by 
the large impact of cosmological surface brightness dimming, 
e.g. by a factor ~150 for z=2.5. This implies that only the high- 
est surface brightness galaxies are spatially resolved with inte- 
gral field spectroscopy or in continuum imaging. It also biases 
the samples where we can detect spatially extended emission 
line gas to galaxies which have intense star formation over many 
tens of kpc^ or extended emissi on line regions excited by p ow- 
erful active galactic nuclei (e.g. jNesvadba et ani2006ll2008h . 

Observations of intensely star- forming galaxies at z~2 with 
bright, extended line emission reveal remarkably broad line 
widths implying important random, and presumably (at least 
partially) turbulent, gas motions. Limited by the ~ 1-few kpc 
resolution of current spectrographs, the observed line widths 
reflect the blended bulk and turbulent gas motions on kpc 
scales. In galaxies where velocity gradients are spatially re- 
solved, both motions have very similar amplitudes, unlike in 
nearby galaxies where line widths are much smaller than cir- 
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cular velocities (Law et^ 20071 2009t iForster Schreiber et alJ 
l2009t lEpinat et al.l l2009l l2010l l2012h This requires either a 
more intense, or perhaps more efficient, source or an addi- 
tional source of turbulent energy at high redshift. The possi- 
ble sources for energizing the turbulence remain the subject 
of significant controversy. Some authors suggest that the tur- 
bulent cascade is primarily initiated by shear and large scale 
gravitational or fluid instabilit ies, perhaps further energized 
by co s mological gas accretion ( Brooks et al.l 20091: Keres et al] 
l2009t iDekel et al.l l2009al: lEknegreen & Burkerd 1201 Ol) . while 
others give a si gnificant role to the energy injection from 
young stars (e.g., Elmegreen & Scalo'^2004t ICoxll2005t iFerraral 
Il993t iNorman & Ferrara 1996) . However, the dichotomy be- 
tween gravitational and stellar sources of energy is largely ar- 
tificial as both no doubt play important roles. Supersonic tur- 
bulence can provide global support to the ISM, but it natu- 
rally produces density enhancements that may collapse locally 
and may subsequently form stars or dissolve. Global star for- 
mation in galaxies appears to be controlled by the balan ce be- 
tween gravity and turbulence (iMac Low & Klessenll2004t) . This 
balance is moderated by the importance of each form of energy 
injection, which on the largest scales depends on the rate and 
scale at which it is injected, and on how rapidly the energy cas- 
cades and dissip ates in various phases, and over what time and 
size scales (^g.- lElmegreen & Scaloll2004l) . None of these vari- 
ables are well constrained in distant galaxies. To some extent, the 
global regulation of star formation mirrors that taking place on 
the smallest scales, where the dissipation of turbulent energy in 
self-gravitating clouds is important for regulating star formation 
wit hin these clouds. 

iLehnert et al.l ( |2009|) found that the velocity dispersions, cr, 
of optical emission-line gas in intensely star-forming galaxies 
at z~2 scale with the star-formation rate per unit area (star- 
formation intensity, Ssfr) as cr °c Sgpp ( see also lSwinbank et all 
l2012tlMenendez-Delmestre et al.ll20i3h . They suggested that at 
a very simple level, such a relationship would be expected if 
the lines were bro adened by the m echanical ener gy injection 
from y oung stars (iDib et al.ll2006h . Subsequentlv. iGreen et alJ 
(1201 Ol) came to similar conclusions from observing strongly star- 
forming galaxies at low redshifts which had similar overall Ha 
luminosit ies as galaxies obser ved at z~2 with integral field spec- 
troscopy (iLehnert et al.l2009|). More recent ly, and through an an- 
alytical approach. lOstriker & She ttvl (1201 lb found that such a re- 



lationship might be the natural outcome of dense (Sgas ~few 100 
Mq pc~^), star-forming gas disks which are regulated through 
the momentum input from star formation. Compared to the other 
possible sources such as gas accretion, star-formation intensity 
in high-z galaxies can be constrained observationally, e.g. by 
measuring Ha luminosities and spatial distributions, although 
heavily extincted regions may have their intensities underesti- 
mated. 

A number of authors subsequently pointed out that this re- 
lationship is_alTO_consistentwith_the Toomre disk stability crite- 
rion (e.g.. lKrumholz & DekelllToiOl) . Although not discussed in 
their work, a Toomre parameter Q~l characteristic of gas disks 
which are marginally unstable against star formation leads to a 
relationship of the type observed, cr oc p ^^^^'', if the s e gala xies 
follow the Schmidt-Kennicutt relation. Burkert et al.l (1201 Ol) de- 
velop this idea further by arguing that disks must stay close to the 
gravitational stability line. They suggest that large scale gravita- 
tional instabilities will generate density and velocity irregulari- 
ties which drive turbulence and heat the gaseous d isk. This pro- 
cess sa turates close to the instability line and thus. lBurkert et al.l 
(1201 ol) argue that gravitational forcing is sufficient to explain the 



highly turbulent gas. JElmegreen & Burkertl(l2010l) find that their 
analytical model of intensely star-forming, fragmenting disks 
can reproduce the observed velocity dispersions in high redshift 
galaxies if Q is set equal to 1, suggesting that whatever the mech- 
anism responsible for driving galaxies towards the line of insta- 
bility, it can reproduce the observations. 

The basic difficulty is that observing velocity gradients and 
line widths consistent with marginally Toomre-stable disks only 
imphes that distant star-forming galaxies are in a state where 
they can form stars intensely, but it does not purport to ex- 
plain the source of the turbulence. Galaxies at z~2 appear 
to maintain their high gas turbulence over long star-formation 
time scales of several 10*^ yrs, and maintaining Q~l over sim- 
ilar timescales is only possible if turbulent dissipation is bal- 
anced by the injection of turbulent energy over similar length 
of time. Since t urbulent energy dissipates rapidly, on only a 
dynamical time (jMac Low! Il999l) . this seems difficult to en- 
sure with gravity alone over t he extended star-formation periods 
(JElmegreen & Burkertll2010l) . Moreover, gravitationally-driven 
processes alone are perhaps insufficient for producing even 
the relatively modest dispersion observe d in the warm neutral 
and ionized medi um of nearby galaxies (JTamburro et al.l 120091 
iLerov et al1l2008l) . Correlations between s tar-formation rate and 
Hi velocity dispersions in nearby galaxies (JTamburro et al.l2009l) 
as well as the decreas e in CO line dispersions with radius in 
some nearby galaxies (IWilson et al.l l201 1) favor an important, 
but perhaps not dominant, contribution from the energy injec- 
tion by supemovae and stellar winds in regulating the turbulence 
of the ISM in nearby galaxies. The same also appears true for 
intensely star-forming relatively nearby galaxie s with properties 
similar to high redshift Lym an break galaxies (JBasu-Zvch et al.l 
l2009l:lGoncalves et alJl2010l) . 

Here we present new observational data on how star forma- 
tion intensities and gas kinematics in high-redshift galaxies are 
related, extending our initial sample from 1 1 galaxies to 53, and 
covering a larger redshift range. Apart from putting our previous 
claims on statistically more robust grounds this also allows us to 
demonstrate the effects of surface-brightness dimming directly. 
We discuss the inevitable biases in all studies of star formation 
in high-redshift galaxies related to surface-brightness dimming 
and beam-smearing from our data, and demonstrate that our re- 
sults are not dominated by t hese effects. With th is larger sam- 
ple we confirm the results of iLehnert et al.l (l2009l) that the ISM 
in the observed high-redshift galaxies has very high pressures 
compared to galaxies at low redshift. We argue that this enables 
star formation in these galaxies to be self -regulated, and demon- 
strate that self-regulation leads very naturally to a Toomre pa- 
rameter, Q~l, generalizing the analytical model first presented 
by lElmegreen & Burkert (2010). We further paint a picture of 
what the optical emission line gas may be showing us about the 
turbulent cascade and its impact on star formation in high red- 
shift galaxies. 

Throughout the paper we adopt a flat Hq =70 km s"' Mpc"^ 
cosmology with Qa - 0.7 and Qm - 0.3. 

2. Source characteristics 

Our sample of 53 galaxies was observed with SINFONI on the 
ESO-VLT. Theirbasic properties of the 53 galaxies are presented 
in Table [1] They span a redshift range of z= 1.3-2. 7 (FiglT). For 
details of the observations of the galaxies in this sample see the 
summaries i n e.g. Lehnert et al. (2009); Forster Schreiber et al.l 
( I2OO6II2OO9I) . All \ yere reduced by us in a similar way as previ- 
ously described bv ILehnert et al.l (l2009l) . 
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Fig. 1. Redshift distribution of our sample of galaxies. The two 
observed redshift ranges of the Ha and [Nii] /16583 lines (z<1.8 
and z>2)are due to atmospheric transmission and the limited 
wavelength coverage of the H- and K-bands. 



The data, which were retrieved from the ESO archives, were 
taken for a variety of programs and the selection of the galaxies 
is rather inhomogeneous. As we are not performing a population 
study, for the purposes of this paper, this is adequate. The key to 
our sample is that all galaxies have sufficiently high Ha surface 
brightness to obtain spatially resolved line maps in a few hours 
integration time (which is typical for SINFONI observations of 
distant galaxies). Because of this, our results are only applica- 
ble to galaxies with high star-formation intensity (star-formation 
rate per unit area) which necessitates a process that is able to 
drive intense star- formation for 100 Myrs or more (lErb et al.l 
|2006|) and my introduce a bias. However, it should be pointed out 
that the characteristics of the ensemble of galaxies studied here 
are similar to those of the general distant star- forming galaxy 
population. Specifically, the galaxies in our sample for which 
the spe cific star formation rate has been estimated (a major- 



ity; see [F orster S chreiber et al.ll2009[ 1201 it iGenzel et al.ll20Tll 



IVergani et al.„2012l and references therein), lie along or near the 
"ma in sequence of star formation" for redshifts between 1 a nd 2 
(e.g. lElbaz et~aDl2007HDaddi et al.ll2007HElbaz et al.ll201 ih . 

The observations reach surface brightness detection limits in 
Ha of ~2xl0"'' erg s"' cm"^ pixel"' (for the 125 mas pixel"' 
scale and averaged over 3x3 pixels; "mas" is milliarcseconds). 
The spectral resolution is FWHM~115 and ~150 km s"' in the 
K and H bands, respectively. The observed surface brightnesses 
range from about 3.6-34 x 10"'** erg cm"^ s"' arcsec"^ and the 
objects have isophotal radii, at their surface brightness detection 
limits, of 1-2 arcsec^, or on average ~7+2 kpc (corresponding to 
an isophotal area of ~ 150+40 kpc^). The point-spread function 
FWHM of the data is ~0''6, which at z=2 represents an area of 
~20 kpc^, so we have generally ~8 spatial resolution elements 
per object. With such low spatial resolution, it is necessary to 
discuss the impact of beam smearing on the characteristics of 
our data (discussed in detail in §|4li. In the present analysis, we 
will formally define the isophotal radius, riso, as riso - (Aiso/tt)'^^, 
where Aiso is the isophotal area, i.e. the total area of all pixels 
above the 3-cr surface brightness limit of the data (Table[T}. 

The total Hff luminosities of all galaxies in our sample are 
above 10"" erg s"' and most are above ^10"*^ erg s"'. Hff lumi- 
nosities in this range imply an average star-formation rate of ~70 
M© yr"'. However, as these are not corrected for extinction our 



estimates are lower limits. For a sub-sample of these galaxies we 
have Hp estimates, and the line ratios of Hff to H/3 suggest rel- 
atively modest ext inction corrections t o the star-formation rates, 
of a factor of few (iLehnert et al.ll2009l) . 

3. Surface brightness limits and selection effects 

Due to cosmological surface brightness dimming, over the red- 
shift range of a factor ~2 in our sample there is about a factor 
of 5 difference in the faintest rest-frame line emission surface 
brightness levels that can be probed, whereas the physical di- 
mensions in the rest frame, i.e., length per arcsec, changes by 
only about 5%. Given this strong impact of cosmological surface 
brightness dimming , it is worth repeating the conclusion from 
iLehnert et al.l ( l2009l) that even the lowest redshift sources in our 
SINFONI sample are extreme compared to galaxies in the local 
volume, where only the most intense starbursts have such high 
Hff surface brightness, and this on considerably smaller physical 
scales only. 

The total range in surface brightness probed by the SINFONI 
observations of our ensemble of galaxies is *20 (Fig. |2]i. For 
purposes of this comparison and others, we constructed 2- 
dimensional histograms which show the frequency of occur- 
rence of various pairs of values in all pixels of our data sets, 
after smearing to the spatial and spectral resolution of the real 
SINFONI observations - for example, surface brightness versus 
projected distance or Ha line width versus star- formation rate 
intensity. Using the frequency of occurrence is a clearer way of 
visualizing various relationships than showing all single pixels 
and emphasizes the continuity of the data. 

To investigate the dependence of the range of surface bright- 
nesses probed by our data and its influence on the observed 
galaxy isophotal size, we divided our galaxies into 4 bins of 
equal galaxy number based on their isophotal areas (Fig. |3)- For 
the 3 bins with the largest average isophotal sizes (>1.4 arcsec^), 
the observed dynamic range in surface brightness goes up to ~20 
and is typically about 15. However, for the apparently smaller 
galaxies, the dynamic range of the observations is only about a 
factor of 3. The relationship between surface brightnesses and 
projected distance is the same for all bins, indicating that their 
individual declines in surface brightness follow the same trend. 
The observations imply two things. First, although the effects of 
beam smearing cannot be ignored, as we will discuss extensively 
in the next sections, it does not dominate the overall distribution 
of the surface brightness (i.e., the distribution does not simply 
represent the point-spread-function, PSF, of the observations). 
Second, it implies that the reason the apparently small galax- 
ies appear small is that the surface brightness detection limits 
of the observations are too high to detect their extended emis- 
sion. As noted above, this is not simply due to the physical scale 
of the galaxies changing across our redshift range, as this effect 
amounts to only ~5%. While we do not show it here, there is a 
crude trend for the galaxies with the highest total Hff fluxes to 
also have the largest isophotal sizes. Again, this suggests that we 
are limited by the overall dynamic range of the data and by the 
impact of surface brightness dimming, but not by differences in 
the intrinsic isophotal sizes of the galaxies (Fig. |4]i. This is an 
important point to bear in mind when reading this paper 

Furthermore, if the trends in surface brightness and total Hff 
luminosity were dominated by cosmological surface brightness 
dimming, which is proportional to ( 1 H-z)"*, we would expect our 
data to show a slope of about 4 in a plot of the logarithmic rela- 
tion between Ha luminosity and isophotal area and indeed, we 
find such a trend (Fig.|4]i. The range in surface brightness detec- 
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Fig. 3. Two-dimensional histograms of the relative frequency of occurrence of normalized Ha surface brightness values versus the 
projected distance for galaxies in 4 equal-size bins in galaxy number selected on their isophotal area (Aiso). The relative frequency 
of occurrence is scaled in the same way as in Fig. |2] From left to right, each panel represents an increasing angular size (i.e., Aiso 
< 1.4 arcsec^, 1.4 arcsec^ >Aiso < 2.5 arcsec", 2.5 arcsec- >Aiso < 3.3 arcsec^ and Aiso ^3.3 arcsec^). The values of riso, which 
is derived from Aiso, for individual galaxies are given in Table [1] The galaxies with the smaller dynamic range (higher surface 
brightness detection limits) in the data are also smaller in isophotal size. 
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Fig. 2. A 2-dimensional histogram of the frequency of occur- 
rence of pixel values of Ha surface brightness versus projected 
physical distance for the entire sample. As shown in the color 
bar, the darkest regions have the highest frequency of occur- 
rence. The surface brightness has not been corrected for cos mo- 
logical dimming. The zero of the projected distance was chosen 
to be the symmetry point in either the Ho- velocity field or of 3-cr 
isophotes in the Ha surface brightness map. 




Fig. 4. Plot of Ha luminosities versus isophotal area for the 
galaxies in our sample (red circles). The size of each circle is 
proportional to the redshift of the galaxy, with larger circles de- 
noting higher redshifts (range in z is 1.3-2.7). The two solid lines 
show the relationship between luminosity and area for the av- 
erage surface brightness detection level of galaxies, 1.1x10"'^ 
erg s"' cm"^ arcsec"" (black for galaxies at z=2.2 and red for 
z=1.3, respectively). The dot-dashed line shows the relationship 
for the minimum surface brightness detection limit of 3.6x10"'** 
erg s"' cm"^ arcsec"^ for all galaxies with za:2.2. See Table [T] 
for the Ha luminosities, isophotal radii and surface brightness 
detection limits. 



tion limits in our data trace out the ranges in size and luminosity 
we observe in our sample of galaxies. In addition, this confirms 
that at the high redshift end of our sample the intrinsic surface 
brightnesses are very high and that they generally decline with 
redshift. 



4. The effect of beam smearing 

In lLehnert et al.l (l2009h . we proposed that the Ha emission line 
velocity dispersions of the high surface brightness spatially- 



resolved galaxies observed with SINFONI (and other near- 
infrared integral field units) could be explained by the mechani- 
cal output of young stellar populations within those galaxies. In 
fact, we framed the argument in the context of an overly simpli- 
fied energy injection model of the form, cr- (eEsFk)'^^, where cr 
is the Ha velocity dispersion, e is the coupling efficiency of the 
mechanical energy output to the interstellar medium and Esfr is 
the star-formation intensity. The energy due to the young star s 
is given by population synthesis models (iLeitherer et al.lll999l) . 
Using scaling relations based on models and observations of 
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nearby galaxies, we were able to explain the trends between dis- 
persion and star-formation intensity with no free parameters (see 
[Lehnert et al. 2009, for details of the scalings used in this anal- 
ysis). 

However, the limited spatial resolution of integral field spec- 
troscopy of z~2 galaxies causes an artificial broadening of the 
emission lines due to large-scale kinematics, which can result in 
a significant radial velocity gradient across each pixel. While we 
argued that this would not have an important impact on our re- 
sults in lLehnert et al.l(l2009h based on some simple modeling, we 
now feel it is important to revisit this issue in a more complete 
and detailed way. 

4.1. Construction of artificial IFU observations 

In order to study the effect of beam smearing, we used 
the simulation of an isolated galaxy already presented in 
IDi Matteo et a l. (200 8J) which was made using a Tree-SPH 



code (see ISemeHn & CombesI l2002t iDi Matteo et"ai1 120091: 



IChilingarian et al.ll2010i for a description). This model has a 



high gas mass fraction (50%) and the galaxy wa s evolved in iso - 
lation, without companions or tidal interactions (lOu et al.ll2011l) . 
We allowed the disk to go unstable against star-formation and it 
developed a "clumpy" morphology which evolved with time (see 
iDi Matteo et al.l l2008l for details of the initial conditions and 
evolution of the simulation). Intrinsically, our simulations show 
low gas velocity dispersions, of-order 10 km s"' (or less), which 
are roughly constant with radius and a rotation speed of ~220 
km s"' which implies a v/crs;20. From this simulation, we pro- 
duced maps of the star-formation intensity, radial velocity field 
and velocity dispersion (Fig. |5]l. For our primary analysis, the 
simulation was viewed at an inclination angle of 45° and scaled 
to have a projected rotational amplitude of 1 10 km s"', a value 
derived from the average intrinsic full width at half maximum 
of the integrated spectra of our sample. This implies a rotation 
speed of about 160 km s"' . 

We then degraded the resolution of the simulations to pro- 
duce synthetic observations similar to what we observed with 
SINFONI. We adopted a single redshift (z=2) for determining 
the scalings in the synthetic observations - the exact redshift 
chosen does not matter much as the physical scale per angu- 
lar projected size (kpc arcsec"') changes little (~5%) over the 
redshift range of our sources. We also took into account the in- 
strumental spectral resolution, the seeing using a FWHM for a 
Gaussian distributed PSF of 0"6 and added an amount of noise 
to the synthetic data consistent with the observations. The artifi- 
cial datacubes created in this way were then analyzed using the 
same procedures as the real SINFONI observations. 

In the following analysis, we will focus on the impact of 
beam smearing on the interpretation of the underlying cause of 
the high dispersion observed in distant galaxies with spatially- 
resolved line emission observations. 

In lLehnert et al. 2009i we proposed that the dispersion is pro- 
portional to the star- formation intensity. For our analysis of the 
effects of beam smearing, we make a comparison between the 
analogs produced using two sets of scaling relations for the ve- 
locity dispersions: (1) the dispersions taken directly from the 
simulation, cr^im, and (2) a function o f the star formation rate 
intensity cr — 140 a/Ssfr km s '- see iLehnert et al.l ( l2009l) for 
details. The resulting Yia flux and kinematics maps made with 
our artificial SINFONI data sets and applying these two different 
types of dispersion show, as expected, that the surface bright- 
ness distribution is now smoother than in the original simula- 



tion (Fig.|6ll. The star-formation intensities - the star-formation 
rate per unit area - are related to the Yia surface brightness by 
a simple scaling (e.g. iKennicutt 1998j). We use this scaling to 
produce the Ha surface brightness distribution which is subse- 
quently beam smeared. 

4.2. Effect of beam smearing on star-formation intensities 

The areas of intrinsically high star formation intensity are 
blended after smearing the data, which lowers their intensities 
(Fig.|6]l. Whether or not such peaks in the star-formation inten- 
sity can be observed in principle depends on their relative dis- 
tribution and peak intensity relative to their surroundings within 
one or two PSFs (cf. Figs. |5]and|6]l. However, Fig. |7] demon- 
strates that, after smearing, the simulated final radial distribution 
of surface brightness is consistent with what we observe. The 
similarity seen does not depend on the underlying velocity dis- 
persion distribution assumed in the construction of the synthetic 
data. 

The star-formation intensities from the simulation decline by 
$25% after beam smearing. While the precise value of this de- 
crease is likely to depend on the initial spatial distribution of 
the surface brightness, it is always going to be small. The total 
star-formation rate, which in the simulations is typically about 
60 Mq yr"', remains essentially unchanged after making an ar- 
tificial SINFONI data set. The largest impact of smearing is on 
the distribution of the star-formation intensities and in losing the 
lowest surface brightness emission because of noise and the rel- 
atively low dynamic range of the data. 

4.3. Effect of beam smearing on velocity dispersions 

Since we have assumed two quite different relationships between 
velocity dispersions and star-formation intensity, it is in the re- 
sulting velocity dispersions that we see the most dramatic dif- 
ferences. In the model with cr^i,^, the dispersion is low (a few to 
10 km s"') and almost constant as a function of radius and star- 
formation intensity (Fig. |8] left panel). In the synthetic obser- 
vations however, the beam smearing by the PSF (FWHM=0"6) 
and the low spectral resolution of the data broadens all lines by 
about ~40 kms"' and the noise in the data increases the scatter 
about this mean rise in the dispersions. Even though we cor- 
rect for the intrinsic resolution of the spectrograph, at relatively 
low S/N we tend to over-estimate the widths of lines that are 
only marginally resolved. Our results suggest that for galaxies 
with roughly constant and low dispersions (~40 km s '), a spec- 
trograph like SINFONI will tend to over-estimate the intrinsic 
widths of the lines but still produce a roughly constant trend of 
dispersion with star-formation intensity. 

For the model where the simulation dispersions are scaled as 
o" = 140 -i/SsFR km s"', there are almost no differences between 
the mean velocity dispersions from the simulation itself and the 
artificial data sets (Fig. [8] right panel). Because of the large range 
of the scaled velocity dispersions, there is little effect of beam 
smearing on the average velocity dispersion. Looking at the 
pixel-by-pixel quantities, we obviously lose the relatively rare 
regions and peaks with the highest velocity dispersions (over 
200 kms"'). However, as suggested by the average and intrin- 
sic dispersions, the decrease in the values (ranges and averages) 
is not large, even after the smoothing to low spatial and relatively 
low spectral resolution. Again, because of the loss of the lowest 
surface brightness regions, we also tend to lose information on 
the areas with the lowest velocity dispersions. However, in spite 
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Fig. 5. Output from the simulation of a star-forming disk galaxy as it would appear at z=2 before it was beam smeared. Shown are 
the distribution of the star-formation rate intensity (left panel). Ha velocity field (center panel) and Ha dispersion map (right panel). 
The galaxy is projected at an inclination angle of 45° relative to face-on. The scale on the axes is in kpc and the kinematic axis is 
horizontal. Despite the different units (kpc vs. arc seconds), we note that the sizes of these panels and those in Fig.|6]are similar, for 
ease of direct comparison. 
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Fig. 6. Results of the simulated data sets of a star-forming disk galaxy at z=2 after beam smearing using different relationships for 
the velocity dispersions. Shown are the distribution of the Ha surface brightness (left), velocity field (middle) and velocity dispersion 
map (right). In the upper row, we have simply taken the velocity dispersions as estimated from the N-body/SPH simulation, while 
for the lower row we assumed a relationship between the velocity dispersion and star-formation intensity of the form cr = 140 a/Xsfr 
km s ' . The circle in the upper right of each plot represents the FWHM of the PSF used to make the artificial data set. The scales 
are in arc seconds. We note that we have scaled the images such that the sizes of these panels and those in Fig. |5] are similar, for 
ease of direct comparison. 
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Fig. 7. Comparison of the relative frequency of occurrence of 
pixel values of normalized per pixel Ha surface brightness dis- 
tributions of the whole sample of galaxies (in shades of red) with 
superimposed an artificial SINFONI data set constructed from an 
N-body/SPH simulation (blue crosses). Both observations and 
simulation follow similar distributions. 



of all of this, the general trend for higher star- formation intensity 
to result in higher dispersion is maintained. 

We also halved the signal-to-noise ratio and used the intrinsic 
rotation speed of the simulated galaxy keeping the inclination 
angle at 45deg. In both cases, the overall trend in the data is 
maintained and the increase in the scatter in the star-formation 
intensity-dispersion plane is not significant. 

4.4. Are the dispersions related to Ha surface briglitness? 

We now compare the observed trend between star-formation in- 
tensity and Ha line velocity dispersion with the models we used 
to make synthetic data sets, including the eff'ects of beam smear- 
ing. Fig. |9] shows that the values of the simulated Ha line dis- 
persion cTsim are too low to explain the large line widths we gen- 
erally observe in our sample of galaxies. On the other hand, the 

1 /2 

simulation using the relationship, cr oc S^pj^, is a good match to 
the overall distribution of the data. We note that more properly, a 
constant term should perhaps be included in the relationship, as 
cr oc Sjpj^ + constant, whic h is both warran ted by observations 
and expected theoretically (iDib et al.ll2006l) . At star formation 
intensities of order 10""^ to 10"^ Mq yr"' kpc"^, the observed 
velocity dispersion of the Ha emission lin e are ~5-20 km s ' 
dFathi et al.ll2007t lErroz -Ferrer et al.ll2012l) . We have chosen to 
ignore this constant in our analysis for two reasons. First, veloc- 
ities this small are below the resolution of the SINFONI data. 
Thus we are unable to probe the characteristics of galaxies with 
low star formation intensities. Second, the star formation inten- 
sities of the galaxies we are observing are several orders of mag- 
nitude greater and this including a small constant term would 
have an insignificant impact on our analysis or conclusions. 

Figure [TO] displays the differences in the Hff velocity as a 
function of the star- formation intensity distribution for the galax- 
ies in two different redshift regimes, z<1.8 and z>2. The key re- 
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Fig. 8. In our analysis of the impact of beam smearing on the ve- 
locity dispersion and star-formation intensity, we adopted two 
relationshipsbetween these two quantities. One taken directly 
from the simulations (left panel), and another where we assumed 
cr oc l^^f^ (right panel). The model velocity dispersions are 
shown in red (these are not lines but many individual points). 
We also show the points (blue) after analyzing the artificial data 
cubes derived from the simulation which include the effects of 
beam smearing. The red + (red x) sign indicates the mean star- 
formation intensity and integrated (original input) velocity dis- 
persion for the intrinsic relation between velocity dispersion and 
star-formation intensity, while the black + (black x) sign has the 
same meaning but now for the beam-smeared relationship. 



suit is that the data are consistent with a simple scaling between 
dispersion and star-formation intensity in both samples - how- 
ever, on average the z~1.5 galaxies have lower dispersions and 
star-formation intensities. Also, the scatter is larger for the lower 
redshift galaxies, which perhaps indicates that a single relation- 
ship of the type cr oc -^Esfr does not explain all of the trends in 
the lower-redshift data as well as it does for the higher redshift 
objects. 

At the low end of the star-formation intensity within indi- 
vidual galaxies and over the ensemble of galaxies with lower 
average star-formation intensities, it is likely that we are over- 
estimating the velocity dispersions, and it could be that the dis- 
persions are becoming roughly constant with radius. We say this, 
because our artificial model-based data cubes suggest that ob- 
served dispersions smaller than ~40-50 kms ' are consistent 
with $10 kms ', a level at which we simply can no longer dif- 
ferentiate between a model where the line widths are propor- 
tional to Xgpjj or one where the line widths are consistent with 
what is observed in the MW and other near by galaxies (see, e.g. 
IWada et al.ll200llDib et al.ll2006tlAgertz et al. 2009). 

In Figures [To] we plotted the simulated data against the nor- 
malized distribution of star formation intensity and velocity dis- 
persion. However, this does not give a sense of the differences 
between galaxies at different redshifts in their distribution of 
star formation intensity and velocity dispersion. The distribu- 
tions of the galaxies in the higher redshi ft regime are consis- 
tent with those of 1 1 galaxies presented in lLehnert et al.l(l2009l) . 
Quantitatively, the ensemble of galaxies with z>1.8 and z<1.8 
have a mean star formation intensity of 0.73+0.09 and 0. 16+0.02 
M0 yr"' kpc"^ respectively. Since both samples contain about 
the same number of galaxies, the width of the distribution of the 
higher redshift sample is significantly larger and the difference 
in wid th between the samples is consistent with iLehnert et al.l 
(I2009h . 

Moreover, correcting individual galaxies in the sample 
for the average extinction (derived from fitting the spectral 
energy distribution with stellar population synthesis models; 
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Fig. 9. Normalized 2-dimensional histogram of the frequency of 
occurrence in the data set of Ho- line widths as function of star 
formation rate intensity. The set of points for each galaxy are 
normalized such that all of the galaxies have the same mean star- 
formation rate intensity (their distribution of the velocity disper- 
sion in each galaxy is unaffected). The mean value was chosen 
to match the mean value of the star-formation intensity in the 
model with a = (eSsFR)'^^ as indicated by the red + sign in 
the right panel of Fig.[|](~0.5 Mq yr"' kpc"^) ^^ j^ Fig.|fl the 
darkest regions represent the highest frequency of occurrence in 
the data set. The green x and blue + symbols represent our simu- 
lated data after beam smearing, assuming that the two quantities 
plotted are related as cr = (eSsFR)'^^ and cr-cTfiim, respectively. 



iForster Schreiber et al.ll20lTl) does not reduce the scatter among 
them. 



4.5. The equivalence of integrated and spatiaiiy resoived 
measurements 

So far, we have studied the spatially resolved relationship be- 
tween the Ha velocity dispersion and t he star-formati o n inten - 
sity as probed by the Ha emission. In iLe Tiran et al.l (12011 ah . 
we argued that the relationship between the integrated Ha lu- 
minosity (or equivalently the total star-formation rate) and sur- 
face brightness-weighted mean dispersion of whole galaxies re- 
veals the same underlying physical mechanism as the spatially 
resolved measurements - namely, that the intense star- formation 
is responsible for the dynamics of the warm ionize d medium as 
probe d through the Ha recombination line (see also lGreen et al.l 

l2oToh . 

To further investigate this relationship, we have tested the 
idea that the intense star-formation in these galaxies can even ex- 
plain the integrated measurements of the total star- formation rate 
and the integrated velocity dispersion. This is important: if we 
can show that also these integrated measurements reveal this un- 
derlying relationship, then a greater wealth of data become avail- 
able for analysis, since slit spectra are more commonly available 
than integral field spectra. 

We first note that we find little difference in velocity dis- 
persion for whole galaxies if we simply sum up the Ha emis- 
sion pixel-by-pixel and estimate their integrated dispersion or if 
we weigh each measurement by surface brightness (or signal-to- 
noise). The difference between both methods is only about 15 
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Fig. 11. Mean Ha line velocity dispersion, o-^e^n (in kms"'), 
versus total Ha luminosity (in erg s"') for different samples of 
local and high redshift galaxies. L ocal galaxies (green triangles) 
are taken from lGreen et al.l (1201 Ol) and the model results for our 
sample are represented by red and blue stars. The blue stars show 
the results where we assign to each pixel in a galaxy an Ha sur- 
face brightness based on the model relationship, cr = 140 a/Ssfr 

km s ' and use the observed local Ha velocity dispersion. We 
then estimate the total Ha luminosity by summing up the pre- 
dicted pixel-by-pixel surface brightnesses, and estimate o-jne^,, by 
weighting the pixel-by-pixel velocity dispersion measurements 
by the Ha surface brightnesses derived from the model. 

km s"', where flux weighting tends to give a lower value and a 
scatter of about 15 km s ' . So compared to the mean integrated 
velocity dispersion of ou r sample, thi s represents a systematic 
offset of only about 10% dLe Tiran et al . 2011a). 

To further test the relevance of a relationship of the form, 
cr=(eSsFR)'^^, where e is the efficiency at which the mechani- 
cal energy from a star-formation intensity Esfr is converted into 
turbulence and bulk flows in the interstellar medium, we applied 
it to the spatially resolved dispersion measurements, to esti- 
mate integrated quantities. Specifically, we used the spatially re- 
solved, pixel-by-pixel Ha line velocity dispersion measurements 
to estimate the Ha surface brightness of each pixel. We then 
used these surface brightness estimates as weights when com- 
bining the individual pixel-by-pixel velocity dispersions, to esti- 
mate the integrated mean dispersion in Ha, (Tmean- We summed 
these estimated surface brightness values to estimate the inte- 
grated Ha fluxes and then converted these to total luminosities 
using our adopted cosmology. The results of this analysis show 
that we can reproduce the actual measurements reasonably well 
considering the crudeness of the model and that we only con- 
sidered one coupling efficiency (Fig. fTTT i. So not only does this 
simple scaling relationship reproduce the spatially resolved mea- 
surements, but it can also reasonably reproduce the integrated 
measurements. Apparently, unlike in many models of the ISM, 
the observed velocity dispersion of the warm ionized medium 
does increase with increasing star-formation intensity and rate 
(iDib et al.,i2006; Ostriker & Shetty 2011i) . 



5. Discussion 

Figure [12] shows that [NII]/Ha ratios are not arbitrary, but scale 
with Ha surface brightness. This has already been noted by 
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Fig. 10. Plots of the normalized frequency of occurrence of Ha velocity dispersion versus star-formation intensity in two redshift 
bins, z<1.8 (left) and z>2 (right) (cf. Fig. |9] for the whole sample). The normalizations are different from that used in Fig. |9] and 
in the left and right panels. For both sub-samples, the average star-formation intensity of each galaxy was shifted to the average 
of its ensemble of galaxies (0.16 M0 yr ' kpc^-for z<1.8 and 0.73 M© yr"' kpc^^and z>2). Note that the two sub-samples tend 
to sample different ranges in both velocity dispersion and star- formation intensity. The green x and blue + symbols represent our 
simulated data after beam smearing, assuming that the two quantities plotted are related as cr = (eZsFR)''" km s"' and cr=crsi„,, 
respectively. The green and blue model points have been shifted to overlap with the data (the location of the modeled points along 
the axis of Zsfr is set by the coupling efficiency of the mechanical energy output of young stars to the ISM). 



iLehnert et al.l (l2009h and is confirmed with our new, larger sam- 
ple. The emission line ratios and Ha surface brightnesses are 
related to the gas pressure, column density, intensity of the radi- 
ation field and the number of line-emitting regions along the line 
of sight within each beam. New, compared to our previous anal- 
ysis of this relationship, is that our expanded sample includes 
galaxies at lower redshift, which in turn will allow us to probe 
pressures, column densities and radiation fields in galaxies with 
less extreme Ho- surface brightnesses. 

Measurements of the [Sii] d oublet ratio in indivi dual galaxies 
in a sub-sample of this sample ( iLehnert et al.ll2009l) and th rough 
a stacking analysis of this sample dLe Tiran et al.ll201 Ibl) . sug- 
gest that the typical densities of the warm ionized gas are 
ne~10s to ~500 cm"-'. To compare these estimates, we ran 
a suite of ph otoionization model s using version 08.01 of the 
Cloudy code (iFerland et al.lll998h to show the relationship be- 
tween the [Nn]/16583/HQ' line ratio and the Ha surface bright- 
ness as a function of the hydrogen internal pressure, column 
density and ionizati on parameters for a single ionized cloud (see 
ILehnert et al.l 120091) . Unfortunately, the interpretation of such 
modeling is ambiguous due to the degeneracy between surface 
brightness, line ratio and the number of individual emission line 
regions per large physical beam of our observations. Given the 
crude sampling of our data, it is likely that the densities are lower 
than predicted by a single cloud model, as would be the ioniza- 
tion parameters. Therefore, for both the higher and lower red- 
shift data sets, the likely densities are certainly <1000 cm""* but 
generally larger than a few times 10 cm"^ and are consistent with 
estimates of the electron density directly measured from the [Sn] 
doublet ratio in individual galaxies and a stacking analysis. 



There is an interesting difference in the regions occupied in 
the [Nii]/16583/Hq' versus Ha surface brightness plane by galax- 
ies at z< 1 .8 and z>2. Taken at face value, it would suggest we are 
probing less dense gas in the lower redshift galaxies. However, 
this is undoubtedly due to probing lower surface brightness gas, 
related to the strong influence of surface brightness dimming, 
which allows us to probe the very diffuse, perhaps even shock 
heated, gas in the lower redshift sample (the regions with rela- 
tively low surface brightness and high ratios of [Nii]/16583/Ha' 
in Fig. [12] likely require a contribution to their ionization from 
shocks). This effect is observational, rather than evolutionary. A 
complete census of both redshift regimes would be necessary to 
show that the average pressure or number of overlapping emis- 
sion line regions (i.e. thickness) in the warm ionized media of 
galaxies is declining with redshift. 

The thermal pressures estimated from the electron densi- 
ties are much higher than in the disk of our Milky Way or 
other nearby normal galaxies. The hydrostatic mid-plane pres- 
sure in the MW is about 10^ "'""^^ K cm"^, while for other local 
spirals it can amount to about 10^ K cm~^, though it is typi- 
cally about 10"*^ K cm^^ (iBhtz & Rosolowskvlllooa and refer- 
ences therein). In the nuclear regions of local starburst galax- 
ies, thermal pressures estimate d from the [Sn] doublet rat io are 
typically about 10*'"^ K cm"^ dLehnert & Heckmanl[T996l) . The 
pressures we estimated from the photo-ionization models are 
consistent with t hose of nearby intensely star-form i ng galaxies 
(•Heckm anetaP 119901: ILehnert & HeckmanI [19961: IWang et al.l 
1998) and measurements of density-sens i tive lin e s at high red- 
shift ([Sn] doublet ratio: iLe Tiran etal.1 [201 Ibl: ILehnert et all 
120091: iNesvadba et al.ll2007l) ! 
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Fig. 12. Ha surface brightness (corrected for cosmological 
surface brightness dimming) versus the logarithm of the 
[Nii]/16583/Hq' line ratio. The galaxies have been separated into 
two redshift bins: z<1.8 (top) and z>2 (bottom). The lines rep- 
resent results for photo-ionization modeling (see text for details) 
for 6 sets of conditions: a range of hydrogen densities, log iih 
(cm"^) =1 (in green), 2 (blue) and 3 (red), and of column densi- 
ties, log Nh (cm"^)=21 (solid lines) and 20 (dashed lines). The 
ionization parameters range from log U=-5 to -1 (increasing 
from left to right along the lines, meaning low ionization pa- 
rameters have relatively high log [Nii]/Hff and low Ha surface 
brightnesses). At low ionization parameters the solid and dashed 
lines for the column densities overlap until the lower column 
density cloud becomes density bounded, at which point the mod- 
els reach very low ratios of [Nii]/16583/Ha' (as shown by the 
sharply dropping dotted lines). For clarity, we only show the log 
Nh (cm"^)=21 and log iih (cm^^) =3 line until log U = -3 and 
for log Nh (cm"^)=21 and log n// (cm"-') =2 until log U = -2. 



By comparing these thermal pressures in the WIM to that 
generated by the radiation pressure of the ionizing radiation 
field, it is possible to constrain th e general distribution of the 
WIM dKrumholz & Matzne3l2009h . Our optical line ratios sug- 



gest that the gas is p redominately photoionized by massive stars 
(iLehnert et al.ll2009l) . To estimate the impact of radiation pres- 
sure on an ionized fluid, it is appropriate to only consider the 
ionizing radiation. We can use the ionization parameter to gauge 
the relative contribution of radiation pressure to the thermal pres- 
sure of the WIM. Simple scaling of the definition of the ion- 
ization parameter suggests Prad/PthenTiai= 2.5xl0"~ U (<hv>/18 
eV)/k(T/10'* K) for T^IQ-* K and <hv>=18 eV, where <hv> is 
the average energy of ionizing photons and k is the Boltzmann 
constant. Since the ionizing parameter is always less than log U 
= -1.6 (and mostly < -3), the radiation pressure is lower than 
the thermal gas pressure throughout the nebulae. Finding that the 
thermal pressure dominates the radiation pressure and that the 
ionization parameters are relatively low, suggests that the WIM 
is widely distributed over the volume of the ISM in our sample 
of galaxies, truly part of the general interstellar medium, and not 
simply near or directly associated with the s tar forming clumps 
or regions (see iKrumholz & Matzneill2009l for a detailed dis- 
cussion of this point). This point is crucial for our subsequent 
analysis as it suggests that the WIM is a tracer of the dynamics 
of the general ISM in distant galaxies and not merely of the Hii 
regions or clumps. 

Other models or relationships have been proposed to ex- 
plain the high dispersion in the warm ionized gas that are di- 
rectly related to the ionization state of the gas. I Wisni o ski et al.l 
(1201 2h and Swinbanket al. (2012) propose that the relationship 
between star formation intensity and velocity dispersion can be 
explained by a relation between the luminosity of a Hydrogen 
recombination line (e.g. H/J) and t he gas velocity dispersion, a s 
observed in local Hii regions (e.g. iTerlevich & Melnicklll98lh . 
The sizes of local Hii regions and the relationship between line 
luminosity and velocity dispersion can be understood simply if 
the ga s within the Hii regions is virialized (ITerlevich & Melnickl 
Il981h . In such a hypothesis, the line widths are mainly due to 
virial motions, implies that the warm ionized medium is in close 
proximity to the star forming clumps (within the gravitational 
sphere of influence). Nearby HII regions, for which such a rela- 
tion is appropriate, generally show ionization pa rameters th at are 
high, log U~ -1 to -2(e.g. Lopez-Hernandez et al.l 120131) . but 
our analysis suggests that distant galaxies generally have much 
lower ionization parameters, implying that their emission line 
gas luminosity is dominated by more widely distributed gas and 
not directly associated with the clumps. 

5.1. Source of the pressure in the WIM 

What are the sources of pressure in these distant galaxies? Is 
the pressure generated by the thermalization of the collective 
mechanical energy output of massive stars through the thermal- 
ization of stellar winds and supernova explosions? Observations 
of nearby galaxies suggest that the low ionization line emit- 
ting region is approximately in pressure equilibrium with the 
X-ray emitting plasma, wit h thermal pressures i n the range 



lO'^ t o 10"' dvne cm~^ (e.g.lHeckman et al.lll990l:lMoran et al.l 



119991: IWestmoquette et aill2007l: IStrickland & Heckmanll2009l) . 
Similar pressures are found in models using typical ini- 
tial conditions expe cted for galaxies exhibiting outflows (e.g. 
ISuchkov et al.lll996l) . If these two phases are in equilibrium then 
we can use a simple model to estimate the pressures in the hot 
X-ray emitting plasma to see if it agrees with our estimates of 
the thermal pressures (densities) in the warm ionized media in 
dis tant galaxies. 



IStrickland & Heckman ( 20091) has shown th at the simple an- 
alytic expressions of lChevaUer & Cleggl(ll985l) for the pressure. 
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temperature and velocity of outflows driven by the intense en- 
ergy injection from young massive stars agree with hydrodynam- 
ical simulations. Adopting an d scaling these simple relations 
from lChevalier & Clegd (1 19851) we estimate the central pressure 
in a starburst region as P^^ 3.5x10"'" (^Mo.26)'^~(fthei-maiE4i.9)'^^ 



R 



SB,350 



Q. dyne cm , where ^ is the mass-loading factor, M| 



0.26 



is the mass return rate in units of 0.26 M0 yr" , ethermai is 
the thermalization efliciency of the stellar winds and super- 
novae, E41 9 is the energy injection rate in units of 7.4x10"*' 
erg s"', RsB,350 is the radius of the starburst in units of 350 
pc and Q is the opening angle of the outflow in sr. The 
adopted scalings for the mass and energy injection rates from 
massive s tars are appropriate for a star formation rate of 1 
M0 yr ' dLeitherer et al.|[l999l) within a radius consistent with 
clump si zes and Ha disk thicknesses in high redshift galaxies 
(iForster Schreiber et al. 201 1 ; Epinat et al. 2012i) . For the radius 
scaling, see I Strickland & Heckmanl(l2009l) . 
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Fig. 13. Ha velocity dispersion, cr as function of stellar mass 
surface density, Sstars (M0 pc"^). The lines represent the Toomre 
criterion Q=l for two rotation speeds (v=100 or 200 km s"') 
and three radii (r=3, 5, 10 kpc), as indicated in the legend. These 
values span th e range of radii and velocity shears observed in 
these galaxies (iForster Schreiber et aDl2006l l2009l) . The steUar 
mass surface densities (red dots) are estimates within one effec- 
tive radius for 6 galaxies in our sample with NICMOS images 
(see text). The red solid line shows the relationship between cen- 
tral velocity dispers ion and stellar mass surface density f or local 
early-type galaxies JShen et al.ll2"00llTavlor et al.ll2010l) . 

Reasonable assumptions for the thermalization efficiency, 
mass loading factor, and mass and energy injection rates, result 
in pressures of 10"*^ to 10"'" dyne cm"^ in our galaxies, similar 
to what is estimated for nearby starburst galaxies. This simple 
modeling suggests that the hot plasma generated by the thermal- 
ization of the mechanical energy output of young stars can gen- 
erate sufficient pressure to approximately match our estimates 
for the thermal pressure in the warm ionized medium. 

One of the assumptions that is needed for the model to be ap- 
propriate is that the energy injection r ate is approximately con- 
stant over an outflow timescale (e.g. Chevalier & Clegglll985t 
ISuchkov et al. 1996 ; Strickla nd & Heckman 2009i). This has 
been shown to be true for local starbursts (e.g. iHeckmanet alj 
Il990t lLehnertetanil999l) . While we do not know the outflow 



timescale of the distant galaxies, since their star-formation in- 
tensities are similar to nearby starburst galaxies (see discus- 
sion in Lehnert et al. 2009), it is logical to assume their out- 
flow timescales would be similar (10s Myr or less). Distant 
galaxies similar to the ones o bserved here (e.g. lErb et al.]|2006l; 
IForster Schreiber et al.l 1201 lb have star formation timescales 
of a: 100s Myrs and are likely much longer than the outflow 
timescales. A reasonable assumption is that the energy injection 
rate is approximately constant in distant gala xies as it appears to 
be in n earby starbursts and thus the model of IChevaher & Clega 
(Il985h is an ap propriate approximation for e stimating the pres- 
sures (see also [Strickland & Heckmanll2009l for a discussion). 



5.2. Mass and energy flow? 

The aforementioned rough equality of the thermal pressure in 
the two ISM phases, the hot 10^-10** K plasma and the 10"* 
K warm ionized medium, is important. It may imply that the 
mechanical energy injected into the ISM of the galaxy couples 
efficiently with the recombining gas and perhaps drive a mass 
and energy fl ow from the WIM to the c old neutral and molec- 
ular medium dGuillard et al.ll2009[ l2010l) . The recombining gas 
is unstable at the pressures we observe and will quickly cool, 
collapse and becom e part of the cold neutral medium, CNM 
(IWolfire et al.l 1 19951) . Since this is expected to happen quickly 
this would imply that t he CNM a cquires the relative kinematics 
of the WIM (Guillard e t al .120121) . If the gas is dusty, this recom- 
bining gas will rapidly turn molecular since the H2 formation 
time scales on dust grains is much shorter than the expected dy- 
namical time_scales_oftherecombining, cooling and collapsing 
clouds (JGuillard et al.ll2009h . The formation of gravitationally 
bound clouds releases gravitational energy as turbulence into the 
molecular gas. 

This picture links the energy input from intense star forma- 
tion to energy dissipation, the formation of molecular gas and 
star formation. Thus not only is this a flow of mass and energy 
but also one of kinetic energy, which will not only contribute 
to the turbulent energy of the molecular gas but also serve to 
balance dissipation. Of course, it is not a simple one-directional 
flow, but a cycle or feedback loops where individual molecu- 
lar clouds may form stars when they become self gravitating, 
or simply be heated and destroyed. Naturally, because all the 
timescales in the cycle are relatively short, <10 Myrs, the ISM 
will be in a dynamical equilibrium. The components of the ISM 
should not be viewed as single, isolated components but rather 
as an intricate balance whose relative mass flow rates are deter- 
mined by a competition between the heating and cooling rates as 
well as the outflow rate due to starburst driven winds and the in- 
fall rate due to the cosmological accretion of gas through cooling 
gas and mergers. 

One of the conditions that must be met in such a picture is 
that the mass deposition rate of the WIM should be more than 
sufficient to fuel star formation, or Msf < estarfoi-mMHn, where 
fstarform is the efficiency of star formation. The recombination 
rate of the Ha emitting gas is given by simple ionization bal- 
ance, Mhii=1.9x10^ L]^q,43 Mq yr"'. The star formation rates 
in these galaxies are ~ 10-300 M0 yr"' and thus the WIM recom- 
bination rates can supply the gas for star formation as long as the 
star formation efficiency is >0. 1 %. For a reasonable estarform only 
a small fraction, a few to about 10 percent, of the recombining 
gas would then be necessary to fuel this cycle. So at least in 
terms of mass flow rates, there is sufficient mass in the warm 
ionized medium to plausibly sustain such a cycle. 
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Of course, for the energy injection to be effective in regu- 
lating star formation, as we have hypothesized, it should drive 
the galaxy towards the line of Toomre disk stability, Q~l. This 
occurs if the star formation can inject sufficient turbulence to 
roughly balance gravitational collapse and hydrostatic pressure. 
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Fig. 14. Ha velocity dispersion, cr, as function of the star- 
formation intensity, Ssfr (Mq yr ' kpc"^). The lines represent 
the Toomre criterion for gas, assuming both Q=l and that the 
Schmidt-Kennicutt relation is appropriate for a variety of rota- 
tion speeds and radii, as in Fig.[T3]and as indicated in the legend 
(see text for details). The underlying display of the normalized 
data (shown in various shades of red) is the same as in Fig. |9] 
(see that figure and text for details). 



While the WIM traces the overall dynamics of the ISM 
(see beginning of § [S]), its contribution to the turbulent pres- 



sure is negligible. It is given by P 



gas,turb 



< Pg 



(T~^ where 



< Pgas > is the volumeweighted mean density (=Pphase/v,phase, 
where pphase is the intrinsic density and /v.phase the volume filling 
factor of the ISM phase). The total turbulent pressure is actually 
the sum of those of all individual ISM phases, but we have sup- 
pressed the summation for clarity. The Wa emitting gas likely 
does not contribute significantly to the total turbulent pressure. 
Assuming simple case B recombination, Tp^lO"* K and equal- 
ity of proton and electron densities, we can crudely calculate the 
emitting volume needed to supply the Ha luminosity compared 
to the volume over which it is observed. The volume filling fac- 
tor of the recombining gas is ffv.wiM ~ Vem/Vn - 4x 10"^ Lj^^, 43 
"e300 ^n 100' whcrc Vem is the emitting volume, V„ 100 is the ac- 
tual ISM volume in units of 100 kpc^ (the typical isophotal area 
is 100 kpc~ and we assumed a disk thickness of 1 kpc), LHa43 
is the Ha luminosity in units of 10'*"' erg s"' and ne.300 is the 
electron density in units of 300 cm""*. The filling factor of the 
recombining gas is ~10"*. 

The most significant contribution to the total turbulent pres- 
sure is likely from the cold neutral and molecular medium. We 
can estimate the density of the dense molecular gas using ob- 
servations of distant galaxies with star formation rates and gen- 
eral properties similar to ours. H2 gas mass surface density 
measurements are avail able for galaxies simil ar to those stud- 
ied here (z~l and 2: iDaddi et al. 2010 ; Arayenaet all 1201 Ot 
Dannerbauer et al.l 120091: iTacconi et al.l 120101: ICombes et al.l 
201 ih which suggest that at z~2 the mean value is about 500 



M|T) pc~^ and that total g as fractions -30-50% (iDaddi et al.l 
120101 : Ir^cconi et aDl2010l) . The volume weighted H2 gas den- 
sity, < riH, >=2gas/H mn, = 10.2 2gas,50o/Hikpc, where Egas.soo is 
the gas mass surface density in units of 500 Mq pc"^ and Hi^pc is 
the thickness of the gas layer in units of 1 kpc. This implies that 
Ptui-b,moi=3 xlO"^ < pgas.io > o'h:,ioo' whcre < ph,,io > is the H2 
volume weighted particle density in units of 10 cm""* and cth,, 100 
is in units of 100 km s"'. These estimates sugge st that the turbu- 
lent p ressure in the molecular gas may be high (ISwinbank et al.l 
l20Tll) . 

For a galaxy to lie near the line of stability, Q=l, the 
pressure due to gravity and that due to turbulence should be 
roughly equal. To estimate the pressure due to gravity we will 
assume hydrostatic equilibrium, which is obviously an over- 
simplification. However, for the present analysis, we are con- 
cerning ourselves with globally averaged estimates, suggesting 
that hydrostatic equilibrium is a reasonable assumption. It is 
impossible to get direct estimates for the hydrostatic gas pres- 
sures for our sample of galaxies as stel lar mass surface densi- 
ties a re available for only a few objects (iForster Schreiber et al.l 
l2011h . As with the estimates of the volume-weighted gas mass 
densities, we will have to rely on results from similar sam- 
ples of distant galaxies, which find stellar mass surface densi- 
ties of 100-2000 Mr:^ pc'^ (eg. iBarden et al.ll2005l: lYuma et al.l 



I2OIII : iForster Schreiber et al.l 1201 11) . The hydrostatic pressure 
scales as Pgas,hydro=7xl0""' Sga.s,50o5;totai,i500 dyne cm"2, where 
2^gas,500 is the gas mass surface density in units of 500 Mq 
pc"- and Stotai.isoo is the stellar mass surface density in units 
of 1500 Mq pc"^. For this estimate of the hydrostatic pressure, 
we have assumed that the velocity dispersion of the gas and 
stars are equal, which is consistent with our general underly- 
ing hypothesis that the velocity dispersions of the warm ionized 
and the cold molecular ISM phases are similar. However, this 
is not a tight constraint as lowering the dispersion in the gas 
would lower the turbulent pressure as well as the hydrostatic 
pressure for a fixed stellar velocity dispersion. The estimated 
ran ges of gas fraction and stellar surface de nsities for z~2 galax- 
ies dPaddi et alJl20IO|: iTacconi et alJl20IOh suggests hydrostatic 
pressures of Pgas,hydro~10"^ to 10"'° dyne cm"^. Thus it appears 
that the hydrostatic pressure and the turbulent pressure in the 
molecular gas are similar if the molecular gas has kinematics 
similar to that of the WIM. 



5.3. Surface densities and tiie Toomre Q parameter 

We have argued that if the turbulence in the molecular gas were 
similar to that of the WIM, then the distant galaxies should 
lie near the line of stability. We proposed a qualitative sce- 
nario where the kinematics of the WIM may be roughly cap- 
tured by the cold molecular gas and thus lead to high turbulence. 
Furthermore, we suggested that t hese intens i ties ca n be under- 
stood in terms of self-regulation (ISilklll997l |200I|) . Such self- 
regulation would reveal itself through mass surface densities in 
the galaxies consistent with Q=I. 

iKrumholz & Dekell (1201 Ol) favor disk instabilities as the 
source of the high velocity dispersion in distant galaxies with 
spatially-resolved ki nematics and line emission (see also, e.g. 
ICacciato et alJl2012h . Their argument is that this is the simplest 
explanation, because it is based on the Toomre instability pa- 
rameter, Q^Kcr/nGl,, where k is the epicyclic frequency (which 
in our subsequent analysis we have taken to be 2'^^v/r, where v is 
the rotation speed and r is the radius) and £ is the mass surface 
density of stars or gas (IWang & Silklll994t iRomeo & Wiegerfl 
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l201lh . They suggest that the high velocity dispersions are a 
natural result of the need to keep the star-formation intense, 
which implies that Q~l. The results given above suggest per- 
haps that the star-format ion itself may also lead to Q~l (see also 
lOstriker & Shettvll20TTb . Combining the Schmidt-Kennicutt re- 
lation with the Toomre instability criterion, Q=l, suggests that 
cr oc Sgpj^. This is similar enough to our suggested relationship 
as to be indistinguishab l e jWisnioski et al. 2012; S winbank et ajj 
I2OI2I) . ISwinbank et alj ( |2012|) show that assuming a marginally 
stable Toomre disk combined with relationships for the largest 
unstable fluctuations in disks and a Jeans relation for the growth 
of unstable modes, can lead to insights into why disks at high 
redshift may have high dispersions, why clumps have high mass 
surface densities, why their masses in disk galaxies decline with 
redshift and why their mass distributions have exponential cut- 
offs. However, such a solution does not explain how the observed 
high dispersions are maintained. This either requires a long dis- 
sipation time or a large source of mechanical energy (iMac Lowl 
Il999t lOstriker & Shett"vll201 ih . What maintains the instabihty 
such that Q~l? We suggest it is the star-formation itself that 
maintains the necessary dispersion in the gas phase, as we will 
now argue. 

The Toomre parameter for the stellar population is 
Q=A:cr/:TGSstais, where Sstais is the mass surface density of stars. 
For the galaxies we have presented here, the velocity shears are 
of-order 100-200 km s"' (Forster Schr eiberet al. 2006, 200j) 
and th e radii are about a few to 10 kpc (iForster Schreiber et alJ 
l2011h . The mass surface density, 2stars=Mstai-s/27r(b/a)r^, where 
Mstars is the total stellar mass, b/a is the minor-to-major axis 
ratio and r^ is the effective radius. Using NICMOS and other 
imaging data of 6 galaxies in our sample, we have made crude 
estimates of stellar mas s surfa ce densities using the results from 
iForster Schreiber et al.l (1201 ll) . We find that the mass surface 
densities within one effective radius are ~100-1500 Mq kpc"^. 
Assuming Q=l, we can relate the velocity dispersion to the mass 
surface density and it appears that the estimated mass surface 
densities are consistent with Q=l (Fig.[T3ll. 

Without considering the effects of beam smearing, this sim- 
ple relation is actually also roughly consistent with our data 
on the emission line properties, as we can demonstrate rather 
simply: if we assum e Q=l for the ga s, and adopt the Schmidt- 
Kennicutt relation dKennicuttl Il998h between star formation 
intensity and gas surface density, we find the relation, cr - 
(7rGQ//c)(EsFR/2.5xlO"'*)'^'^. For a variety of rotation speeds and 
disk radii, we find general consistency with our data (Fig. [T4l l. 
Of course these estimates are crude and could be fortuitous 
but they do lend additional weight to our overall arguments. 
We have not used the generalized (effective) Toomre crite- 
rion, 1/O.ff = 1/Op . + 1/O.t.n,-. (Wang & Silk 1994; Jog 199^ 
iRomeo & Wiegerg|201 ll) because the crudeness of our analysis 
does not warrant this additional complication. However, since 
this is a relatively small change in the adopted Q, it will make 
little difference to our conclusions which is that the galaxies gen- 
erally he near the fine of stability. 

5.4. Toomre-stability of star-forming high-z disks 

Can we make a self-consistent model whereby intense star for- 
mation is responsible for driving high turbulence and push- 
ing the ISM in distant gala xies towards the line of stability? 
lElmegreen & Burkerj (120101) developed an analytical model to 
follow the evolution of gravitational instabilities and star for- 
mation in disk galaxies that are built from rapid gas accretion 
at high, constant rates. The energy of the accretion is balanced 



by the dissipation rate and loss of energy due to star forma- 
tion. The accretion rate is defined by accreting half of a fiducial 
mass scaling in one dynam ical (orbital) time (about 280 Myrs, 
lElmegreen & BurkertI 120101) . Star formation only serves as an 
energy sink in their formalism, to account for gas lost to the 
ISM because it is bound in stars, but they do not include the 
energy return from young stars into the ISM in a self-consistent 
way. In their model, turbulent velocities drop to very low values 
well before stellar-to-gas mass ratios reach observed values. To 
overcome this, they force turbulent velocities, cr, to have values 
consistent with Q= 1 . They find that ultimately the disk reaches 
an equilibrium where Maccretion = Ms,a,-foi-mation- 
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Fig. 15. Six normalized quantities (see legend) as a function of 
time, where t/t' = 90 Myrs or about 1/3 of a dynamical time of 
the potential. To limit the total mass so that it does not become 
too large compared to the observations, the accretion ended at 
t'=6 (about 540 Myrs or two dynamical times). The top panel 
is for a relatively high assumed efficiency of coupling of the 
mechanical energy output of the stellar population to the ISM 
(12%), whereas the bottom panel is for a relatively low coupling 
efficiency (5%). The scalings for these normalized quantities are 
provided in the text (see .Elmegreen & Burkerti2010.) . 

We will now extend their model to include the effects of the 
energy release of star formation. This is necessary if we want 
to investigate how star formation is regulated in galaxies whose 
evolution is not only driven by accretion and dissipation of tur- 
bulent energy, but where star formation contributes significantly 
to the gas energetics, as indicated by our observational finding 
that star formation produces pressures similar to the mid-plane 
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pressure of the disk. In our model, accretion, dissipation and the 
energy released by star formation together set Q. Including star- 
formation fe edback in this way is one sig nificant difference to 
the model of lElmegreen & BurkertI (l20I0l) . The other is that we 
allow for a cut-off in the gas accretion at some point in the evolu- 
tion. Some cut-off or severe reduction in the initial high rates of 
gas supply must have occurred to produce galaxies with masses 
similar to those generally observed at high redshift. This is not 
merely an academic issue either It is often argued that the star 
formation rate and the gas accretion must be equal, and at the 
high star formation rates observed in z~2 galaxies, the star for- 
mation time is a few 100 Myr to 1 Gyr To account for the en- 
ergy from star formation, we have modified equation (11) from 
lElmegreen & BurkertI (|2010|) as follows: 



dE' 



3/2p,I/2. 



,„,(e'm:j'/^-s'e7m;,,+6sfEsfS'm;,,^''^e'"^(1) 



where E', M^jj,|^, M' and S' are the scaled energy, disk 
mass, gas mass and star-formation rate, respectively (see 
[Elmegreen & Burkert 2010, for details). The first 3 terms are 
straight from lElmegreen & BurkertI (l2010l) . The first term is con- 
stant, because for simplicity we assume a constant potential, and 
it is unity, because the energy is normalized relative to the to- 
tal accretion energy, which is extracted from the potential. The 
second term quantifies the turbulent dissipation. The dissipation 
rate is Ecr/H, where E is the total energy, cr is the turbulent 
velocity and H is the scale height of the disk, which is deter- 
mined by a competition between the turbulent energy and the 
gravitational force due to the total mass surface density, i.e., 
H^o-^/nGEdisk- The third term is the energy that is lost to the 
ISM when stars form. Adding a time dependence to the poten- 
tial is quite straight-forward and we have conducted both a time 
independe nt and time dependent analy sis. To these first three 
terms from 'Elm egreen & BurkertI (|20 1 Ob we have added another 
term to incorporate the energy injection rate from star formation, 
espEspS'IVI' ^^^E''''^, where esp is the coupling efficiency of the 
mechanical energy to the interstellar medium of the galaxy and 
Esp is the rate at which energy is ejected per unit of star for- 
mation. The whole equation is normalized such that the accre- 
tion energy is one. In addition, other than the parameters already 
mentioned, other key parameters in determining the evolution of 
this simple model are the star formation efficiency, the gas accre- 
tion rate, the fraction of the accretion energy that couples to the 
kinematics of the ISM, the fraction of mechanical energy pro- 
duced by young massive stars that cycles back into the ISM and 
the dissipati on timescale. For thes e param eters, we adopted the 
values from Elmegree n & BurkertI (|20 1 Ol) . although we allowed 
for a star formation efficiency of 3% in the models instead of 1 % 
they adopted. The fraction of energy produced by massive stars 
that dynamically coupled to the ISM was a free parameter, which 
we would like to constrain, and the dissipation timescale was de- 
termined by the parameters of the model, as just discussed. 

The scalings within this model are quite important. Using 
the estimates of the m echanical energy per unit star formation 
(iLeitherer et al.ll999h and a rotation speed of 220 km s ' (which 
sets the potential energy of the galaxy from which the accretion 
energy is extracted), we estimate that the mechanical energy in- 
jected into the ISM by stars is several 100 times greater than 
that generated by accretion. Also the time scale for star forma- 
tion is much shorter than that for accretion. Thus star formation 
can play an important role in injecting energy into the gas and 
perhaps even pushing disks towards marginal stability against 
star formation (Q~l). The scalings for the original model (taken 



from lElmegreen & BurkertI I20T0I) are t/t' = 90 Myrs (or about 
one-third of an orbital time), S/S' = 160 Mq yr"', Mgas/Mg^^ = 

1.4 xlO'" M© and crlcr'-lQ km s"'. The scaled Toomre param- 
eter, Q' is related to the actual Toomre parameters by a factor of 
^1/2 to 1, depending on the efficiency at which accretion energy 
and the mechanical output from stars is converted into turbu- 
lence. 

The results of this model are shown in Fig. [15] When the 
coupling efficiency of the mechanical energy from the young 
stellar population is high, around 5-15%, Q~l through most 
of the evolution and because of this, the star-formation rate is 
roughly constant. The energy injection into the ISM from the ac- 
cretion and star formation is high and therefore, dispersions are 
also high throughout much of the evolution. With the adopted 
scalings, the dispersions are about 50-150 km s"', and star- 
formation rates are about 50-200 Mq yr"'. All of these values 
are within the range observed in our sample. Star formation in- 
creases the Toomre parameter to Q~l throughout the evolution 
and stabilizes disk star formation. At higher coupling efficien- 
cies, the gas dispe rsion can be even higher but the star forma- 
tion is suppressed. lElmegreen & BurkertI (l20I0b set Q=l when 
Q in their model dropped below one. However, since we allow 
the energy injection from star formation to evolve with star for- 
m ation, we find a different tim e evolution. Unlike in the model 
of lElmegreen & BurkertI (|20 1 Ob . the star- formation rate does not 
reach the accretion rate until well after the accretion has stopped 
and the gas fraction has dropped below 50% of the total (gas plus 
stellar) mass. 

When the coupling with the mechanical energy output of 
the stellar population is weak, the disk is no longer stabilized 
near Q~I and Q becomes v ery small (similar to the results of 
lElmegreen & BurkertI 120101) . Consequently, star formation in- 
creases rapidly and exceeds the accretion rate. This high star- 
formation rate occurs again when the gas fraction drops be- 
low about 50%. In this model, star- formation rates are ~200 
Mq yr"'. This is high but reasonable compared to the star- 
formation rates in our galaxies. The values in Table [1] are gen- 
erally lower than this, but they have not been corrected for 
extinct ion, which is hke ly to increase them by factors of a 
few (jLehnert et al.l l2009l) . We note that during the early, gas- 
rich phase of evolution in these models, the Jeans mass is 
relatively high which implies that the disks would be very 
clumpy, especially if the mechanical energy from stars cou- 
ples very efficientl y to the gas (because su ch galaxies stay 
relativ ely gas rich; Elmegreen & Scalol 20041: Elmegreen et aLJ 
2009allbt iBournaud et al.l 120071: iBournaud & Elmegreenl 120091: 



Ceverinoetal.ll20I0h . 



The basic point illustrated by this simple model is that star 
formation can play a key role in regulating large-scale instabil- 
ities within the disks of distant galaxies. Previous studies have 
come to the same conclusion in a somewhat differ ent context 
(e.g.. lsiia[T997l I2OOIL I2OOI ISilk & Norman! l2009l) . What we 
have added to the analysis is the use of the star formation to 
drive the line widths (which are a combination of turbulent and 
bulk motions). In addition, given a sufficiently high coupling ef- 
ficiency, young stars can control the subsequent star formation 
in distant galaxies. There does not need to be correspondence 
between the instantaneous cosmological gas accretion rate and 
the star-formation rate. In fact, quite the opposite. Smooth con- 
stant accretion is probably unrealistic (e.g. lDanovich et al.l20I2l) 
and what we have shown here is that star formation, even if ini- 
tiated by the accretion event, does not have to be related closely 
in time with accretion. The amount of energy available from the 
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supemovae and stellar winds is much larger in comparison and 
very early in the evolution begins to control the star-formation 
history. 

5.5. Is the coupling unrealistically efficient? 

There are of course the important questions of whether or not 
coupling efficiencies as high as those that are needed to ex- 
plain the large line widths are realistic and if the energy output 
from the young stellar population through SNe and stellar winds 
can drive the necessary turbulen ce. Some modeling studies sug- 
gest t his is not the case (iDekel et alj|2009tt^ Krumholz & Dekell 
|20i3), while other suggest that it mav be (iJoung et al.ll2009h . 

Observationally, there are fe w constraints on what the value 
of the coupling factor might be. 'Krumholz & Dekel (1201 Ol) crit- 
icized the analysis of Lehneit et al. (2009), arguing that the re- 
quired coupling efficiency between the mechanical energy from 
star-formation and the ambient ISM is too high. Our analysis 
in the previous section suggested that coupling efficiencies of 
5-15% are necessary to explain the approximate characteristics 
of distant galaxies, especially the relatively high velocity disper- 
sions. But even these estimates are very uncertain in that we do 
not really know the dissipation time scales and that changes in 
the dissipation rate will modify the coupling factor. They must 
be short, approximately one turnov er time scale of turb ulence 
(a few to tens of Myrs), as argued in lLehnert et al.l(l2009l) and as 
suggested by the model discussed in the previous secti on and the 
simulations and physical argume nts presented in (e .g. iMac Low! 
[199 9: Elmegreen & Scalo 2004; Joung et al.ir2009h . The model 
suggests that the dissipation time scale decreases with time and 
increases with coupling efficiency. Since the dissipation rate is 
proportional to the ratio of the velocity dispersion and the scale 
height of the disk, the disk must be thick throughout much of its 
evolution. This agrees with our finding that the gas has a rela- 
tively low ionization (and weak radiation pressure) so it must be 
widely distributed in the ISM. 

While we call this the "coupling factor", there are various 
factors that contribute to this coupling. To unwind these, we 
will use the well studied nucleus of M82 as the analog for our 
distant galaxies, given that the star formation intensity of M82 
on a several 100 pc scale is similar to that of the galaxies pre- 
sented here. First, the mass ejection from the stellar population 
must be efficiently thermalized. In a study of the nuclear X- 
ray emission properties, which probes gas up to temperatures 
of ab out lO^^K and directly probes the "piston" driving out- 
flows, [Strickland & HeckmanI (l2009l) find that on scales of kpc, 
the thermalization efficiency is 0.3-1.0. The entrainment factor- 
the amount of ambient ISM that is swept-up, accelerated, and 
mixed into the hot hig h pressure plasma - is also i mportant. 
Again, the modeling of IStrickland & HeckmanI ( |2009|) suggests 
relatively low entrainment rates of 1-4 times the mass loss of 
the stellar population in the region surrounding the nucleus of 
M82. In addition, in the optical emission li ne gas of M82, large 
line w idths and high pressures are observed (IWestmoquette et aH 
l2009l) . This suggests that a fraction of the energy from the in- 
tense star formation is dissipated in the WIM. 

There are several ways in which the coupling between the 
thermalized mechanical energy from young stars may be en- 
hanced. Star formation which is distributed increases the cou- 
pling factor betwe en the mechanical energy and the cold /warm 
ambient medium dFragile et al.ll2004t ICooper et aT1l2009l) . This 
enhancement results from the fact that with distributed regions of 
star formation, say in a clumpy disk, overlapping bubbles of hot 
plasma begin to interact and cause the formation of dense rims 



of gas in this interaction zone (iTenorio-Tagle et aDl2003ir2006l) . 
This also has the effect of reducing the opening angles through 
which the hot plasma can flow (JTenorio-Tagle et al. 2003) which 
would thereby increase the central pressures in the star form- 
ing regions, further enhancing its influence on the surrounding 
ISM. Our photoionization and simple dynamical modeling im- 
plies that our galaxies have a geometrically "thick" (~1 kpc) 
warm ionized medium of moderate densities and low ionization 
parameters. Such a thick medium would enhance the interaction 
between the mechanical luminosity of the young stars and the 
ambient ISM. The strong interaction zones, the reduced opening 
angle of the outflows, and a thick wann ionized medium would 
have the effect of smothering the outflow and increasing the ef- 
ficiency at which the mechanical energy couples to the ambient 
ISM, and also likely increases the cooling of the hot plasma gen- 
erated by the collective actio n of SNe and stellar winds in each 
star fo rming clump or region (IWiinsch et al.l2008l iHopkins et al.l 
l20I2h . There is evidence that the optical er nission line gas iri 
distant galaxies may be geometrically thick (lEpinat et al.ll2012l 
e.g.). In addition, these interacting and overlapping regions of 
strong mechanical energy input may also drive gas to collapse 
and form m olecular clouds which may be sites of subsequent sta r 
formation (iTenorio-Tagle et al.l 120051; iHennebelle et al.l |2008|) . 
thus aiding in completing one part of the gas cycle we are ad- 
vocating. 



Furthermore, we have up to now only considered the influ- 
ence of the mechanical energy output of the massive stars but 
not their radiation field (except in the WIM itself). Depending 
on the star formation history and duration, a significant frac- 
tion of up to many times the momentum imparted by the 
mechanical output of a sta rburst could be imparted on dusty 
gas b y radiation pressure (iLeitherer et al.l Il999[ iMurrav et al.l 
l2005l) . Moreover, if the optical depth to IR photons is high, 
the radiation pressure m ay dominate the total momentum out- 
put dHopkins et al.ll201ll) . However, it remains difficult to judge 
the importance of radiation pressure in the overall pressure 
(cfjMurrav et al. 2005; Silk & Nusse^|2010l lOstriker & Shettvl 
1201 ll : iHopkins et alJl201 ll) . For our discussion, the important 
point is that radiation pressure is more likely to be effective in 
imparting momentum to the dense, dusty molecular gas whereas 
the mechanical energy output will tend to flow around, engulf, 
and destroy the molecular clouds. Thus radiation pressure may 
be very effective in driving the large scale kinematics of the 
denser gas and thus regulate star formation globally. 



The coupling efficiency needed to drive the turbulence is 
dependent on the energy dissipation timescale of the gas. In 
the range of coupling efficiencies that lead to reasonable results 
compared to observations, the energy dissipation timescale is in 
the range of a few to tens of Myrs or in other words, shorter than 
a dynamical time or the duration of the intense star formation, 
and approximately as long as one turnover times cale of turbu- 
lence for a 1 kpc thick Ha emitting layer (see iLehnert et al.l 
|2009|) . Thus in such a model, the cycle proposed is predicted 
to be long-lived through repeated cycles of a turbulent cascade 
driven by the energy input of massive stars. Since the rotational 
lines of H2 are one of the dominant coolants of warm molecu- 
lar gas, and trace the dissipation of supersonic turbulence (e.g. 
Guillard et al. 2010), observations of distant galaxies with MIRI 
on the JWST in the infrared will be crucial in constraining the 
coupling efficiency and testing the scenario we have proposed. 
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5.6. Summary 

Through our analysis of near-infrared integral field spectro- 
scopic observations we found that the optical emission line gas 
has a high thermal pressure compared to nearby galaxies, and 
ionization parameters that range from those of giant Hn regions 
to those of diffuse low density gas away from star formation re- 
gions, similar to what is observed in riearby normal and star- 
burs t galaxies (e.g. Lehnert & Heckman 1996;'Wan2 et al.'l998^, 
Le Tiran et al. 2011b; Kewlev & Dopita 2002 ; Pellegrini et al.i 
201 ih . The broad optical emission lines observed (cr]^^<50-250 



km s"') are not mainly due to beam smearing but show that 
lines are intrinsically broader in regions of more intense star 
formation. This suggests that the warm ionized gas has a high 
level of turbulence and bulk motions, and par ticipates in out- 
flows (iLe Tiran et al.ll201 IbHOenzel et al.ll201 ih fliat are largely 
driven by the intense star formation. The high level of turbulence 
and the influence of intense star form ation is consistent with the 
clumpy nature of distan t galaxies (e.g. lElmegreen & Scalol2004i 



; or aistantg; 
an i2009allM) 



lEknegreen et 

The total pressure of the warm ionized medium appears to 
be dominated by the thermal pressure as its average volume- 
weighted density is very low. The relatively low ionization pa- 
rameters of the emission line gas suggest the thermal pressure in 
the WIM is due to the mechanical energy output of massive stars 
and not radiation pressure of the ionizing field. 

The evidence we have presented favors the high turbu- 
lence being driven by mechanical output of young stars, 
which is the energy source that regulates the gas to be near 
the line of disk instability (Q~l). The mechanical energy 
from the intense star-formatio n is likely efficiently thermalized 
dStrickland & HeckmanI l2009h . creating large regions of high 
turbulence. The intense radiation field from young stars may also 
impart momentum to the gas, depending on the gas column den- 
sities and opacity to IR photons. The importance of the radiation 
pressure is that it may couple to the dense dusty gas more effec- 
tively than the mechanical energy output, and thus regulate star 
formation globally. 

Both the mechanical energy (through the thermalized hot 
plasma generated by stellar winds and supemovae) and the ra- 
diation field of young stars generate turbulence, bulk motions 
and high turbulent pressures in the gas which are similar to 
the thermal pressure of the WIM and the hydrostatic pressure. 
However, for the turbulent pressure to be similar to the hydro- 
static pressure, the high velocity dispersions observed in the 
WIM have to be (at least partially) transferred to the dense 
molecular gas, which likely has a much larger mass and aver- 
age volume weighted density than the recombining gas. Such 
transfer may happen because at high pressures the recombining 
gas is unst able and will quickly cool and collapse to a cold neu- 
tral phase (IWolfire et al.llI995l) . and become molecular if dusty 
dGuillard et al.ll20()9l) . Since this collapse and cooling time scale 
is much shorter than the dynamical time sc ale, the kinematics o f 
the WIM and the CNM should be similar dGuillard et al.ll2012l) . 
This leads to the formation of gravitationally bound clouds 
which have dissipated their gravitational energy as turbulence 
in the molecular gas. Some of these clouds will form stars, and 
some may be destroyed by the intense energy output of mas- 
sive stars. The ISM will quickly reach a dynamical equilibrium, 
due to the relatively short timescales for gas cooling, and the 
formation of molecular gas and gravitationally bound clouds. 
The short timescale compared to those of star formation indi- 
cates there is sufficient time for repeated cycles, allowing for the 
build-up of the stellar population in our distant galaxies. 



To explain the observed star formation intensity - line width 
distribution, we made a simple model for the energy flow in a 
galaxy undergoing both intense gas accretion and star forma- 
tion. This model suggests that only a few to perhaps 15% of the 
injected energy from stars is required to drive the high level of 
turbulence we observe. This relatively small fraction of the total 
energy injection is reasonable. 

If the high dispersion observed in the optical emission line 
gas is due to turbulence and if the denser phases of the ISM share 
such high turbulence, then the gas and stellar mass surface den- 
sities in z~2 galaxies are consistent with Q~ 1 and star formation 
is a plausible hypothesis for driving the gas towards the line of 
stability. If so, this would be an avenue for self-regulated star 
formation in distant galaxies. More observations, especially of 
the warm molecular gas in distant intensely star forming galax- 
ies are needed to gauge whether or not this picture is correct. 
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Table 1. Properties of the High-z Galaxies 
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[NII]i6583 


1.3822+0.0009 
1.3824+0.0009 


85+ 13 
83+28 


1.48+0.22 
0.54+0.18 


13+3 


6.0 


6.7 


DEEP2-32015501 


10800 


Ha 

[NII]^6583 


1.3919+0.0007 
1.3921+0.0006 


58+5 
32+ 12 


2.47+0.18 
0.38+0.14 


23+2 


3.8 


9.5 


DEEP2-32021317 


5400 


Ha 

[NII]i6583 


1.3815+0.0017 
1.3822+0.0028 


198+ 20 
348+ 39 


2.92+ 0.28 
3.34+ 0.37 


26+3 


4.3 


8.3 


DEEP2-32021394 


5400 


Ha 

[NII]^6583 


1.3735+0.0009 
1.3753+0.0013 


92+9 
156+ 34 


2.78+ 0.25 
1.29+0.28 


25+3 


4.0 


9.5 


DEEP2-32029850 


6600 


Ha 

[NII]i6583 


1.3952+0.0009 
1.3960+0.0010 


81+22 
108+ 36 


0.13+0.03 
0.09+ 0.03 


1+ 1 


3.6 


2.5 


DEEP2-32037003 


5400 


Ha 

[NII]^6583 


1.3985+0.0014 
1.3987+0.0010 


157+7 
103+ 15 


5.50+ 0.24 
1.40+0.19 


52+3 


3.6 


9.6 


ECDFS- 10525 


5400 


Ha 

[NII]i6583 


2.0261+0.0029 
2.0271+0.0014 


278+ 27 
118+20 


6.90+ 0.67 
2.59+ 0.43 


163+ 16 


11.5 


7.8 


ECDFS-2896 


3600 


Ha 

[NII].l6583 


2.3060+0.0017 
2.3050+0.0019 


137+ 17 
155+ 29 


1.55+0.19 
1.14+0.21 


50+7 


16.0 


3.5 


F257 


21600 


Ha 

[NII]i6583 


2.0264+0.0014 
2.0255+0.0011 


122+ 11 
89+29 


1.92+0.16 
0.38+0.12 


45+4 


2.8 


9.3 


K20-ID4 


9900 


Ha 

[NII]^6583 


1.5934+0.0008 
1.5937+0.0008 


62+4 
64+7 


2.85+0.14 
0.79+ 0.08 


37+2 


1.7 


8.3 


K20-ID5 


12000 


Ha 

[NII]i6583 


2.2240+0.0022 
2.2240+0.0029 


195+ 10 
260+ 23 


2.97+0.14 
1.83+0.16 


88+5 


5.6 


6.1 


K20-ID7 


25200 


Ha 

[NII]^6583 


2.2233+0.0017 
2.2243+0.0013 


148+7 
111+ 14 


4.38+0.18 
0.87+0.11 


130+6 


2.9 


9.4 


K20-ID8 


15000 


Ha 

[NII]i6583 


2.2230+0.0011 
2.2241+0.0015 


89+5 
123+ 17 


2.20+0.11 
0.74+0.10 


65+4 


3.7 


7.5 


Q1623-BX376 


17400 


Ha 

[NII]^6583 


2.4080+0.0011 
2.4085+0.0007 


83+6 
32+ 11 


0.36+ 0.02 
0.05+ 0.02 


12+ 1 


4.6 


3.5 


Q1623-BX447 


13800 


Ha 

[NII]i6583 


2.1471+0.0014 
2.1498+0.0005 


118+ 11 


1.07+0.10 
<0.35'' 


29+3 


4.2 


5.8 


Q1623-BX455 


12000 


Ha 

[NII]/16583 


2.4065+0.0014 
2.4080+0.0015 


107+8 
124+ 34 


1.41+0.09 
0.35+ 0.09 


50+4 


10.7 


4.2 


Q1623-BX543 


8400 


Ha 

[NII]i6583 


2.5204+0.0017 
2.5160+0.0011 


137+6 

72+21 


5.94+ 0.26 
0.68+0.19 


239+ 11 


12.5 


6.1 


Q1623-BX599 


5400 


Ha 

[NII]^6583 


2.3305+0.0019 
2.3302+0.0011 


160+7 
85+ 17 


5.99+ 0.24 
0.76+0.15 


199+9 


8.2 


6.2 


Q1623-BX663 


33300 


Ha 

[NII]i6583 


2.4253+0.0026 
2.4265+0.0038 


219+ 12 

325+ 58 


1.96+0.11 
0.75+0.13 


72+5 


4.3 


6.2 
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Table 1. continued. 



Object 


Int 


Line 


z 


cr 


Flux 


SFR 


SB limit 


Tiso 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Q2343-BX389 


14400 


Ho- 

[NII]i6583 


2.1716+0.0020 
2.1727±0.0014 


180+7 

122+ 27 


4.59+0.17 
0.69+0.15 


128+5 


3.4 


8.5 


Q2343-BX502 


9600 


Ho- 

[NII]^6583 


2.1550±0.0008 
2.1571±0.0010 


52+3 


1.82+0.09 
0.23'' 


50+3 


2.8 


7.7 


Q2343-BX528 


34200 


Ho- 

[NII]i6583 


2.2684±0.0015 
2.2691±0.0010 


128+6 
72+ 10 


2.26+0.10 
0.41+ 0.05 


70+4 


4.2 


7.5 


Q2343-BX6I0 


21600 


Ha 

[NII]^6583 


2.2098±0.0018 
2.2105±0.0013 


163+7 
106+6 


7.91+0.32 
2.49+0.13 


231+ 10 


3.7 


9.3 


Q2346-BX482 


31200 


Hff 

[NII]i6583 


2.2562+0.0015 
2.2566+0.0010 


131+6 
70+ 10 


4.09+0.17 
0.48+ 0.06 


126+6 


3.4 


9.0 


SA12-5241 


7200 


Ha 

[NII]^6583 


1.3617+0.0007 
1.3619+0.0006 


53+3 
37+4 


1.46+0.07 
0.62+ 0.06 


13+ 1 


4.3 


6.0 


SA12-6192 


10800 


Ha 

[NII]i6583 


1.5038+0.0008 
1.5038+0.0009 


62+5 
84+ 12 


0.69+ 0.05 
0.40+ 0.06 


7+ 1 


1.8 


5.7 


SA12-6339 


19200 


Ha 

[NII]^6583 


2.2963+0.0013 
2.2950+0.0016 


102+5 
129+ 26 


3.33+0.15 
0.61+0.12 


107+5 


4.3 


6.0 


SA12-8768 


10800 


Ha 

[NII]i6583 


2.1875+0.0010 
2.1875+0.0011 


73+4 
85+21 


2.31+0.12 
0.44+0.10 


66+4 


5.5 


6.6 


SA15-5365 


9000 


Ha 

[NII]^6583 


1.5340+0.0011 
1.5323+0.0010 


117+ 12 
100+ 40 


1.10+0.10 
0.19+0.07 


13+2 


2.1 


6.8 


SINS-4751 


3600 


Ha 

[NII]i6583 


2.2651+0.0011 
2.2651+0.0010 


78+4 
69+ 17 


3.12+0.16 
0.48+0.11 


97+5 


6.2 


5.8 


SSA-22a-MD041 


21600 


Ha 

[NII]^6583 


2.1632+0.0014 
2.1627+0.0016 


122+6 
142+ 39 


3.54+0.16 
0.50+0.13 


98+5 


2.9 


8.3 


VVDS-020147106 


7200 


Ha 

[NII]^6583 


1.5182+0.0008 
1.5183+0.0007 


76+4 
40+ 13 


5.37+ 0.27 
0.31+0.10 


62+4 


4.1 


7.6 


VVDS-020261328 


3600 


Ha 

[NII]^6583 


1.5276+0.0008 
1.5279+0.0006 


71+4 
39+ 13 


1.00+0.06 
0.12+0.04 


11+ 1 


7.9 


4.3 


VVDS-220014252 


7200 


Ha 

[NII]i6583 


1.3090+0.0011 
1.3090+0.0011 


130+6 
128+ 35 


7.84+ 0.35 
1.02+0.28 


63+3 


5.6 


8.5 


VVDS-2200 15726 


10800 


Ha 

[NII]^6583 


1.2920+0.0012 
1.2919+0.0010 


140+7 
111+ 17 


6.63+ 0.30 
1.09+0.16 


51+3 


5.2 


7.7 


VVDS-220544103 


7200 


Ha 

[NII]i6583 


1.3955+0.0010 
1.3960+0.0012 


103+5 
141+23 


7.77+ 0.35 
1.55+0.25 


73+4 


5.6 


8.7 


VVDS-220584167 


7200 


Ha 

[NII]^6583 


1.4643+0.0012 
1.4645+0.0010 


125+6 
98+ 10 


9.37+ 0.42 
1.94+0.18 


99+5 


2.6 


11.5 


VVDS-220596913 


10800 


Ha 

[NII]i6583 


1.2649+0.0010 
1.2654+0.0011 


107+5 
125+31 


4.11+0.18 
0.70+0.17 


30+2 


3.3 


8.5 


ZC-1101592 


3600 


Ha 

[NII]i6583 


1.4034+0.0013 
1.4035+0.0006'' 


150+ 16 

7+2'' 


1.07+0.11 
0.24+ 0.07'' 


10+2 


10.6 


4.1 


20-782941* 


13800 


Ha 

[NII]/16583 


2.1813+0.0016 
2.1814+0.0010 


144+6 
77+ 10 


3.14+0.13 
0.74+ 0.09 


89+4 


3.1 


4.1 


20-782941" 


7200 


Ha 

[NII]i6583 


2.1811+0.0015 
2.1812+0.0010 


133+6 
76+ 19 


4.97+ 0.21 
0.71+0.17 


140+6 


6.2 


7.2 



Table 1. Column (1) - Object designation. " denotes a 250 mas pixel scale plus adoptive optics, and * a 100 mas pixel scale plus adaptive optics. No 
indication denotes that the data were taken without the benefit of AO, and at 250 mas pixel"' ; Column (2) - Integration time in seconds; Column 
(3) - Line identifications; Column (4) - Redshifts of the Ha and [NII]/16583 lines in the integrated spectrum. By integrated spectrum, we mean 
the sum of the fluxes from each pixel with an Ha line signal-to-noise ratio >3 in a specific data cube (see column 8 for the 3cr surface brightness 
detection limits). For each object, the sum for the [NII]/16583 emission line was determined over the same aperture as for Ha; Column (5) - Ha 
line velocity dispersion of the integrated spectrum of each galaxy, corrected for instrumental resolution and in units of km s"' ; Column (6) - Ha 
and [NII]/i6583 line fluxes of the integrated spectrum in units of 10"'* erg s"' cm"-. "^The [Nil] line is not detected and the value given is the 3-(T 
upper limit, for a line width equal to that of the Ho- line, and ''The estimated line characteristics are influ enced significant ly by a nearby night sky 
line; Column (7) - Star Formation Rate, in Mq yr"', estimated from the total Ha luminosity following JKenrricutlll 1 9981 ) . These values have not 
been corrected for extinction, which is likely to increase them by a factor of a few (see text); Column (8) - Ha surface brightness detection limit 
in units of 10"" erg s"' cm"^, at a signal-to-noise level of S/Na;3. These are average values over a 3 x 3 pixel area; Column (9) - Isophotal radius, 
riso, defined as Ajso = n r^^^. The isophotal area, Ais^, is defined as the projected area on the sky above a signal-to-noise level of 3 in the Ha data. 
riso is in kpc. 
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